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Abstract 

An electrochemical Ni-PU battery model has been expanded to include thermal effects. A 
thermal energy conservation equation was derived from first principles. An 
electrochemical and thermal coupled model was created by the addition of this equation 
to an existing multiphase, electrochemical model. Charging at various rates was 
investigated and the results validated against experimental data. Reaction currents, 
pressure changes, temperature profiles, and concentration variations within the cell are 
predicted numerically and compared with available data and theory. 

Introduction 

The nickel hydrogen battery provides an important form of energy storage. This battery 
has applications in both the aerospace and the electric vehicle industries. Nickel 
hydrogen batteries are used in these fields due to their ability to provide long cycle 
lifetimes, high specific energy, high power density, and high tolerance to overcharge. 

Demands by industry have made mathematical modeling of these batteries crucial, in 
particular, thermal modeling has become essential. Temperature variations within the 
battery influence the performance, life, and reliability of the power source. At low 
temperatures there are capacity losses due to high internal resistance within the battery 
and at high temperatures the capacity is compromised by fast self-discharge. Use of the 
active material must be balanced in order to maximize battery life. To do this a uniform 
temperature profile is necessary. Modeling provides a way to monitor the temperature 
profiles, and a way to test improvements without running lengthy and expensive 
experiments. Another very important reason thermal modeling is essential is safety. Side 
reactions occurring within the battery often result in temperature increases, which could 
lead to thermal runaway. Modeling would allow the runaway to be predicted. 

The cell electrochemical behavior is described by electrochemical reactions that are 
affected by temperature-dependent electrochemical and transport properties. This 
interrelation of the electrochemical reactions and temperature-dependent properties 
makes it necessary to couple the electrochemical and thermal models of the NLH 2 cell 
into a single comprehensive model. Incorporating the thermal energy conservation 
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equation, in which the local heat generation rate is determined by the electrochemical 
processes, into the model achieves the coupling. 

A one-dimensional, thermal and electrochemical model has been developed by DeVidt et. 
A1 [2]. The transport phenomena of the gaseous species and electrolyte within the cell 
were described using concentrated solution theory and the volume averaging technique. 

A pseudo-second dimension was used to model the proton diffusion inside the active 
material. Three reactions in the nickel electrode were considered. The main reaction was 
the oxidation/reduction of the nickel species into the solid active material, and the side 
reactions were oxygen evolution/reduction and oxidation of dissolved hydrogen. The 
present work provides a two-dimensional, fully coupled, thermal and electrochemical 
model. It also makes use of the volume averaging technique and concentrated solution 
theory [8]. Unlike the model developed by DeVidt et. al. [2], the present work employs a 
multiphase reaction scheme. At the nickel electrode two main reactions are considered. 
Oxidation/reduction reactions of the nickel species in both the beta and gamma phases 
are considered. The side reactions, oxygen evolution/reduction and oxidation of dissolved 
hydrogen, are also included. The gamma phase is less stable than the beta phase and 
occurs mainly during overcharge [6]. Adding both the beta and gamma phase reactions to 
the model allows for more accurate prediction of the discharge/charge behavior as was 
demonstrated by Timmerman [7]. 

DeVidt’s work [2] considers the cell to be isothermal spatially. An adiabatic case and an 
isothermal case were investigated. The focus of this study is to create an axisymmetric, 
multidimensional, thermal and electrochemical coupled model. The multidimensional 
nature of the present model is necessary because temperature gradients in the Ni-Fh cell 
are present mainly in the radial direction due to heat dissipation through the battery wall, 
while the electrochemical and ion transport processes occur in the axial direction. 
Temperature gradients within the cell will be determined by this model due to their 
importance in predicting cell performance, making improvements, and for safety 
purposes. The energy equation of the two-dimensional model will include a heat 
generation term that incorporates the heat effects due to electrochemical reactions, joule 
heating, and phase changes of non-electrochemical nature. Also, as mentioned earlier, 
DeVidt’s model [2] considered two extreme boundary conditions, an adiabatic case and 
an isothermal case. This model allows for a general convective boundary condition. 

Thermal-Electrochemical Model 

A single Ni-H2 cell consists of an aluminum current collector, a nickel electrode, a 
separator, a hydrogen electrode (which has a micro-porous Teflon membrane backing 
permeable only to gas), a nickel current collector, and a gas diffusion screen, as shown 
schematically in Fig. 1 . A concentrated KOH aqueous solution serves as the electrolyte. 
The Ni-H 2 battery is composed of a stack of these disk-shaped cells housed within a 
pressure vessel [3]. The cells are center aligned within the vessel. There is an outer gap 
between the wall of the pressure vessel and the circumference of the cell stack to ensure 
electric insulation between the cell stack and the pressure vessel and to provide a path for 
hydrogen gas to diffuse from the external reservoir into the cell stack. It is important to 
note the cell is axisymmetric, therefore although only the radial and axial directions will 
be examined the results will be the same as those found for a three-dimensional case. 
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Figure 1. Schematic Representation of a Ni-FF Cell 


Assumptions 

The present work is based on the following assumptions. This problem is considered to 
be axisymmetric. The convection effect between the outer edge of the cell and the 
pressure vessel wall, important when the temperature difference between the two is large, 
is ignored for simplicity. As a result, in this space the species transport is by diffusion and 
migration only and the heat transfer is by conduction only. The nickel electrode consists 
of composite cylindrical needles with a substrate inside (see Fig. 1). The active material 
film is assumed to be a mixture of multiple phases with a constant diffusion coefficient 
for the protons. Furthermore, the a and y phases are treated as a single phase [7]. The 
electrode porosities are considered to be constant. There is a continuous gas-pore network 
throughout the cell with a uniform constant volume fraction. The gas, a mixture of 
hydrogen, oxygen, and water vapor, is ideal. Its composition varies with time but does 
not change with position [1], The solid phase is assumed to be completely wetted by the 
electrolyte film. Interfacial chemical equilibrium exists in the liquid phase for all species 
other than the dissolved oxygen that has a small value of mass diffusivity in the liquid 
electrolyte. Electrical equilibrium exists in the liquid phase due to the large value of ionic 
conductivity of the electrolyte. Finally it can be said that the combined thickness of the 
current collector and the Ni electrode is much smaller than the diameter of the cell, 
therefore <f> SU b is considered to be constant along x but to vary with r. 

Model Equations 

At the nickel electrode the following main reactions occur: 

(5 phase reaction: 
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mOOH p + H 2 0 + e~ <» Ni(OH) 2 „ + OH~ 
y phase reaction: 

NiOOH r + H 2 0 + ne~ <=> Ni(OH) 2 a + OH ' 
y phase production: 


(la) 


(lb) 


NiOOHp + ne~ + K + -> NiOOH r + O/T 


(lc) 


a phase reconversion: 

M(0//) 2 . a -> M(Oi50 v 


(Id) 


Oxygen evolution/reduction and the oxidation of dissolved hydrogen also occur at the 
nickel electrode and are considered to be side reactions: 


-0 2 +H 2 0 + 2e- <=>20H- 
2 2 2 


2H 2 0 + 2e~ <— H 2 + 20//' 


( 2 ) 

(3) 


Two reactions occur at the hydrogen electrode, the main reaction: 


H 2 +20H- o2// 2 0 + 2e‘ 


(4) 


and the side reaction is: 


2 OH' <--0 2 +H 2 0 + 2e 
2 2 2 


(5) 


The general Butler-Volmer equation is used to determine the rates at which these 
reactions occur. The equations are written with respect to a specified reference state, 
namely a temperature of 10°C and an electrolyte concentration of 7. 1 d-3 mol/cm 3 . 
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The subscripts in the above equations; la, lb, lc, 2, 3, 4, 5 refer to the reactions labeled 
with the respective equation number. 


The overpotential of a given reaction j, qj, is given by [4]: 

( ( \[ dU dU Kf \\ 

'l J = K-t-\U J , re/ + (T-Tj—- 


cT dT 


(13) 


The temperature dependency of the overpotential is considered in the above equation. 
The reference overpotentials for the nickel electrode are calculated using the following 
equations [7]: 
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and the reference overpotentials for the oxygen and hydrogen reactions are [1]: 


RT ( \ 

U>,ef = U 4j ref = - 0.9263 — ln(/7"; ) (16) 

u 2 „, = l/ w = 0.3027 + — \r{p % ) (1 7) 


k in equation 14 is an intercalation constant used to describe the interaction between 
species in a solid solution. 


The conservation of OH species in the liquid phase is represented by the following 
equation [7]: 
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where t°_ is the transference number of OH' with respect to the velocity of the solvent. 
j OH is the total transfer current from all electrochemical reactions that generate or 
consume OH at the electrode/electrolyte interface: 


f. OH .OH 
J b +J g 


in the nickel electrode 


j 0H = S 


0 


in the separator 


< aH(T n 4 + Tn5) in the hydrogen electrode 


(19) 


where 


j OH p- aNi(l - X y )[T nla + 0 n2 + t n3)p] (20) 


j = a Ni X y [I n 1 b + (l n2 + T n3) y ] 

Similarly, j 2 and j 2 are the transfer currents associated with the generation or 
consumption of hydrogen and oxygen respectively: 
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in the nickel electrode 
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( 21 ) 
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Conservation of H 2 species in the liquid phase is given by [7]: 

d ( £eC " 2 ^ = V(£# Vcf 2 )+ 2- j H > + J" g 


dt 
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O 2 species conservation in the liquid phase is [7]: 
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The interfacial mass transfer rate of the hydrogen or oxygen from the liquid to the gas 
phase is expressed as: 


Hi 


RT C g ' C ej 


i = H 2 and 02 


J = K» 

eg V 

Conservation of species in the gas phase is given as [7]: 




dt 


Conservation of species in the solid phase [7]: 
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where 


aNi(l -Xg)T n la for5 = b 

aNiXglnlb for5 = g 

Conservation of charge in the liquid phase is represented by [7]: 

V ■(k eff V<t > e )+ v(kf V In c OH )+ j OH =0 
Conservation of charge in the solid phase is given [7]: 
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The energy conservation equation is written in the following form [4]: 
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The third term of the energy equation has three components. The first part: 
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(32) 


pKZ',fc+ n>r 

represents the heat effect produced by the electrochemical reactions [4]. This component 
includes both the reversible heat characterized by the Peltier coefficient: 

dU , 

n ,=T — J - (33) 

J dt 

and the irreversible heat related to the surface overpotential. 

The second component of the heat generation term is associated with phase changes of 
the non-electrochemical nature, e.g. evaporation or condensation of water within the cell 
[4]: 


£A/,r = AVv + ^<"‘ 0 - / ."' < ’ 

The third component of the heat generation term: 


(34) 


7T +kfvinc 0l < Vt,)dV (35) 

* y a 

represents the joule heating in both the solid and liquid phases [4], 

Electrochemical and transport properties that are temperature dependent are expected to 
follow the Arrhenius equation [4]: 


O = <t> 0 exp 
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where E act is the activation energy corresponding to a property O. 

The top and bottom of the cell were considered to be adiabatic. The thermal boundary 
condition applied to the inner edge of the cell and the outer wall of the pressure vessel 
was the following: 

~ k % = h(<T ~ ) ( 3? ) 
where k is the thermal conductivity of the material, and h is the heat transfer coefficient. 

Results and Discussion 

Cell performance under various conditions was simulated using the model described 
above. Physical properties and electrochemical parameters were provided from 
experimental work conducted at JPL and are listed in Table I. Several parameters were 
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not available, but found through trial and error to produce the best overall agreement with 
experimental charge data (these values are indicated with an asterisk). 


Table 1: Parameters 


Parameter 

Values 

References 

Kinetics 

io,la 

3.0d-04 A/cm 2 

* 

io,lb 

1 .0d-04 A/cm 2 


io,ic 

2.0d-55 A/cm 2 


io,2 

3.0d-10 A/cm 2 

* 

io,3 

2.5d-16 A/cm 2 

* 

io,4 

6.5d-05 A/cm 2 

* 

io,5 

1 ,0d-07 A/cm 2 

* 

aia 

5.0d-01 A/cm 2 


aib 

5.0d-01 A/cm 2 


au 

5.0d-01 A/cm 2 



1 ,0d+00 A/cm 2 


a 3 

5.0d-01 A/cm 2 


a4 

5.0d-01 A/cm 2 


0.5 

1 ,0d+00 A/cm 2 


P hydrogen, ref 

52.0 Atm 

Ref 1 

P oxygen, ref 

5.0e-4 Atm 

Ref 1 

U, a 

4.2d-01 V 


u Ib 

3.35d-l V 


U, c 

4.4d-01 V 


u 2 

3.0d-01 V 


U 3 

-9.26d-l V 


u 4 

-9.26d-l V 


u 5 

3.0d-01 V 


Ela 

5.0d+04 J/mol 


Eib 

3.0d+04 J/mol 


Eic 

0.0d+00 J/mol 


e 2 

8.0d+04 J/mol 


e 3 

0.0d+00 J/mol 


e 4 

2.0d+04 J/mol 


e 5 

1 .2d+05 J/mol 


Thermal 

dudtu 

-2.1 Id- 3 V/K 


dudtib 

-1.51 d-3 V/K 


dudtic 

0.0d+00 V/K 


dudt 2 

-1.68d-3 V/K 


dudt 3 

-8.36d-3 V/K 


dudt* 

-8.36d-3 V/K 


dudts 

-1.68d-3 V/K 


initial temperature 

15°C 


reference temperature 

15°C 


environmental temperature 

0°C 
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heat transfer coefficient 

5 W/m 2 K 

Cell 


electrode pairs 

17 

inner pressure vessel radius 

4.435 cm 

outer electrode radius 

4.400 cm 

inner electrode radius 

0.350 cm 

A1 current collector thickness 

0.0015 cm 

Ni electrode thickness 

0.453 cm 

Separator thickness 

0.210 cm 

Pt H 2 electrode thickness 

0.400 cm 

Ni current collector thickness 

0.0010 cm 

Gas screen thickness 

0.0610 cm 

Porosity of substrate 

0.85 

Porosity of Ni electrode 

0.41 

Porosity of Pt H 2 electrode 

0.85 

Porosity of separator 

0.85 

Specific interfacial area of 


substrate 

2000 cm 2 /cm 3 

Specific interfacial area of 


nickel electrode 

3919 cm 2 /cm 3 

Specific interfacial area of 


Pt-H 2 electrode 

2000 cm 2 /cm 3 

H 2 reference concentration 

1.461d-3 mol/cm 

O 2 reference concentration 

1.529d-3 mol/cm 

Reference concentration 

7. 1 d-3 mol/cm 3 


Experimental data for charging at various rates and a constant temperature of 1 0°C were 
used to validate the present model. Figure 2 compares the experimental cell potential 
curves with those predicted. The general shape compared well. The quick rise during 
initial charging was captured, and the second plateau was reproduced. At lower rates 
there is an under-prediction of the potential and as the rate increases there is a slight 
overprediction. Because a single set of parameters was kept for all simulations the 
observed discrepancies were expected. 

Figure 3 compares cell pressure variations measured experimentally with the present 
model predictions. It plots the cell’s pressure against time. Excellent agreement is shown 
for all three charging rates. The change in cell pressure is approximately linear until the 
overcharge period begins. At that point the pressure rise becomes strongly non-linear. 
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Fig. 2: Comparison of predicted and measured charge curves for a Ni-Hh cell at 
different charge rates and constant 10°C. The lines are simulation results, and the 
symbols denote experimental data. 



Figure 3. Comparison of experimental and simulated pressure variations for 
charging at various rates at 10°C. 
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The following results obtained from the computer simulation are intended to further 
reveal the internal processes of the Ni-H 2 cell during charge. 

Figure 4 plots the reaction current distributions within the nickel electrode when the cell 
is being charged at a rate of C/1 and a constant temperature of 10°C to 50% charge input. 
As discussed earlier there are both primary reactions and secondary reactions that occur 
in the nickel electrode. Theoretically, prior to reaching 100% charge input the primary 
reactions are the principle reactions. Beyond 1 00% charge input the secondary reactions 
take over as the major reactions. Figure 4 clearly shows that for 50% state of charge the 
primary reactions, la and lb, are dominant. It is also important to note that these 
reactions are stronger at the electrode/electrolyte interface than at the back of the 
electrode. Inspection of the chemical reactions, la and lb, reveals that OH' is a reactant 
during charging. A large supply of OH' explains why the reaction is more intense at this 
location. 



Figure 4: Reaction current distribution inside the nickel electrode for a cell charged 
at C/1 at constant 10°C to 50 % charge input. 

Figure 5 also shows the reaction current distribution within the nickel electrode when the 
cell is being charged at C/1 and kept at 10°C. It is charged to 150% charge input in this 
case. As expected, the simulation shows that the primary reactions diminish to 
approximately zero and the secondary reactions develop once the battery has been 
charged to 100% charge input. The oxygen and hydrogen reactions, 2 and 3, become the 
dominant reactions once the 100% point is passed. Also, from Fig. 5 it can be noted that 
the secondary reaction rates are larger at the electrode/electrolyte for the same reason as 
stated previously. 
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Figure 5: Reaction current distribution within the nickel electrode of a cell charged 
at C/1 at constant 10°C to 150 % charge input. 

Figure 6 plots the reaction currents in the nickel electrode as a function of the charge 
input relative to the cell capacity. The primary reactions are seen to dominate until 100% 
charge input is reached, beyond that point the secondary reactions take over. Also to be 
noted is that the total current prior to and after 1 00% charge input remains the same. 
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Figure 6. Reaction currents within a cell charged at C/1 and constant 10°C as a 
function of percent charge input. 
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Figure 7 plots the partial pressures of the hydrogen and oxygen within a cell charged at 
C/1 and 10°C. As expected, the increase of pressure due to the hydrogen reaction is 
dominant. The primary reaction at the hydrogen electrode is the main contributor prior to 
1 00 % charge input, and beyond 1 00 % the secondary hydrogen reaction at the nickel 
electrode is adding to the pressure increase. Also at 100 % charge input and beyond the 
secondary oxygen reaction at the nickel electrode begins, producing a rise in pressure and 
contributing to the overall increase in pressure. 



Time [hr] 


Figure 7: Partial pressures within a cell charged at C/1 and 10°C to 150% of the 
cell capacity as a function of time. 

Figure 8 explores the thermal effects within the cell, showing the temperature profile 
within the cell with respect to the cell radius. As mentioned earlier the cell is shaped like 
a disk with a hole at its center. The cell is hottest near the inner edge of the hole and the 
temperature decreases as the outer diameter of the cell is approached, this is expected 
based on the boundary conditions. A temperature change of 3°C is seen. A temperature 
difference of this magnitude influences the kinetic rates of the reactions by as much as 
10 %. 
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Figure 8: Average temperature of the cell charged at C/1 to 150% of its capacity as 
a function of the radius with h = 5 W/m 2 K. 



Figure 9: Average temperature of a cell charged at C/1 to 150% of its capacity as a 
function of time, when h = 5 W/m 2 K and Tj„j t j a i = 15°C. 
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Figure 9 is a graph of the change in average cell temperature with time for a cell charged 
at C/1. This plot shows an initial temperature drop during the early stages and then an 
increase in temperature during the later times. At the end of the charging period the 
battery is hotter than it was initially. The primary reactions dominate when the charging 
time is less than one hour. It is known that the primary reactions are endothermic [5], 
therefore explaining the temperature drop during the beginning stages of charging. The 
hydrogen reaction at the hydrogen electrode is exothermic as is the hydrogen side 
reaction at the nickel electrode. The hydrogen reactions coupled with resistive heating 
cause the temperature rise as charging progresses [2]. 

Experimental temperature measurements were made at the pressure vessel seam. This 
data was then plotted versus the charging time, as shown in Figure 10. An initial drop in 
temperature is followed by a strong increase in temperature. Qualitatively Figure 9, the 
simulation results for the average temperature of the cell, demonstrates a trend very 
similar to the temperature profile found through external measurements as seen in Figure 
10 . 



Figure 10: Experimental pressure vessel temperatures plotted as a function of the 
charging time for a cell charged at C/1. 

Figure 1 1 addresses two important points. It plots the cell proton concentration contour 
for a cell charged at C/1 to 50 % of the cell capacity. It demonstrates the necessity of 
having two-dimensionality in the model. The plot shows that the proton concentration 
varies in both the axial and radial directions. It also shows the effect of temperature on 
the cell. The proton depletion is greater at the inner diameter of the cell, this is related to 
the higher temperatures at this location, clearly demonstrating the thermal, 
electrochemical link. Also, as was noted earlier, there is larger proton concentration 
depletion along the electrolyte interface due to the plentitude of OH'. 
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Figure 11: Proton concentration in a cell charged at C/1 to 50 % the cell capacity. 


Conclusions: 

A two-dimensional, fully coupled, thermal and electrochemical Ni-H2 model has been 
developed. Comparison with experimental data validates the model’s results. This model 
complements experimentation by allowing the user to understand occurrences within the 
cell for various conditions and for many different charging rates. Examples of the 
program capabilities have been given. Work is underway to expand this model to the 
entire cell stack, and the thermal interactions between the cells that arise due to this 
change will be addressed and investigated. 
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List of Symbols 

a specific interfacial area, cm 2 /cm 3 

c' volume-averaging concentration of species i over a phase, mol/cm 3 

C p specific heat, J/g/K 

D' diffusion coefficient of species i in a phase, cm 2 /s 

E ac t activation energy, J/mol 

F Faraday’s constant, 96,487 C/m ol 

h heat transfer coefficient, W/m 2 K 

i n j area-averaged transfer current density of reaction j. A/cm 2 

ipj.ref exchange current density of reaction j at reference conditions, A/cm 2 
j' reaction current density due to production or consumption of species i, 

A/cm 2 

k thermal conductivity, W/mK 

q’” heat generation term 

R universal gas constant, 8.3143 J/mol K 

R S b electrical resistance from solid/substrate interface to the bulk of solid, 

Q cm 2 
t time, s 

t°. transference number of OH' with respect to the solvent velocity 

T absolute temperature of the cell system, K 

Uj, re f open-circuit potential for reaction j at reference conditions measured with 
respect to a Hg/HgO reference electrode, V 
Vj volume of a phase, cm 3 

Greek Symbols 

a a j, a C j anodic and cathodic transfer coefficients for reaction j 
Ah enthalpy change 

£i volume fraction of a phase in the reference volume fraction 

r|j surface overpotential of electrode reaction, j, V 

<f> general symbol representing a property at a temperature T 

X. thermal conductivity, W/m 2 K 

n Peltier coefficient 

<|> potential in a phase, V 

p density of a species I, g/cm 3 

a thermal conductivity W/m 2 K 
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